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ABSTRACT 

According to some models, there may be a significant population of radio pulsars in 
the Galactic center. In principle, a beam from one of these pulsars could pass close to the 
supermassive black hole (SMBH) at the center, be deflected, and be detected by Earth 
telescopes. Such a configuration would be an unprecedented probe of the properties 
of spacetime in the moderate- to strong-field regime of the SMBH. We present here 
background on the problem, and approximations for the probability of detection of such 
beams. We conclude that detection is marginally probable with current telescopes, but 
that telescopes that will be operating in the near future, with an appropriate multiyear 
observational program, will have a good chance of detecting a beam deflected by the 
SMBH. 
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1. Introduction 

Near-infrared observations of stars near our Galaxy's central supermassive black hole (SMBH) 
have revealed a larger number of young, massive stars than can be explained by t ypical star forma - 
tion models (IGhez et al.ll2005l : lEisenhauer et al.ll2005l ). This "paradox of youth" (jGhez et al.ll2005l ) 
has pointed to the development of a possible continuing top-heavy initial mass fun ction (IMF) in 
the region near the central SMBH (jManess et al.1 120071 : iNayakshin &; Sunyaevll2005l ). A top-heavy 
IMF near Sgr A* would imply the existence of a large number of neutron stars in close proximity to 
the cent ral SMBH. Curren t estimates suggest that there could be > 10 4 neutron stars within 1 pc of 
Sgr A* (|Muno et alJbood ) and possibly ~ 1000 pulsars within ~ 0.01 pc of Sgr A* (|pfahl fe Loeb 
20041 ). 



Theoretical predic tions have also suggested a cusp of 



parsec of our Galaxy (jFreitag et al.l 120061 : iHopman &: A lcxandc 



20, 000 n eutron stars in the central 



X-ra y observations have 



been consistent with this number of neutron stars (jDeegan fc Navakshinl 120071 ) but appear to have 
ruled out the presenc e of a larger n u mber of neutron stars (e.g., 40,000). For this paper, we will 
take the results from iPfahl Loebl (|2004l ) and will assume the existence of 1000 pulsars within 
0.01 pc of Sgr A*. We will also follow Pfahl & Loeb in assuming that n(r), the density of pulsars 
as a function of distance from the Galactic center, falls off as r~ 3 / 2 , so that 



n 



— x 10 b pc 



-3 



(r/lpc)- 3/2 



(1) 
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With this assumption it will turn out that the quantity of greatest importance to our predictions 
is the number of pulsars at distances from the central SMBH between 0.001 pc and 0.01 pc. 

In this paper we will consider the possibility that an appropriate program to monitor pulsar 
beams from Sgr A* would detect a beam that is strongly deflected by the central SMBH. Such a 
system would be of great interest, as precision timing of the radio pulses from such a system would 
measure the properties of the spacetime thr ough which they pro pagate. Preliminary work in the case 



of Schwarzschild black holes, presented in I Wang et al.l (|2009bl ) and IWang et al.l (|2009al ) (hereafter 



Paper I and Paper II), has revealed a rich structure and multiplicity of pulses observed in such 
geometries; subsequent work will look at how pulse timing can be used to measure properties such 
as the mass and spin of the SMBH, and measure or constrain deviations in the higher multipoles 
of the spacetime from the predictions of general relativity. In the present paper we will focus on 
determining the likelihood of observing a pulsar in such a configuration. To do this we will first 
calculate the probability that a single pulsar, in orbit around the central SMBH, emits a signal in 
such a way that it is strongly deflected by the central SMBH, and reaches the Earth. From this 
we then infer the probability that the signal from one of the assumed number of pulsars is strongly 
deflected, reaches the Earth, and is detectable by radio telescopes. 

The paper continues in Sec. [2] with a discussion of the model for the pulsar-SMBH system, and 
with assumptions about pulsar characteristics and telescope sensitivities. The heart of the paper is 
the method of computation of probabilities in Sec. [3j Numerical estimates of probability, based on 
this approach, are given in Sec. HI and considerations for an observing program are given in Sec. EJ 
In Sec. [6] we conclude that an observing program, even with current radio telescopes, would have 
some chance of detecting strongly-bent pulsar emissions, while later generations of telescopes will 
significantly increase the likelihood of observing these fascinating systems. 



2. Background, model, and assumptions 



As a simplification in our probability estimates we take the SMBH in Sgr A* to be a Schwarzschild 
hole. It is essentially certain, of course, that the SMBH is rotating, but the angular momentum J 
is currently thought to be only about half of its maximum possible value of GM 2 /c where M is the 
mass of the SMBH (IMelia et al.ll200ll : iGenzel et al.l l2003). Typically, the astrophysical properties 
of Kerr holes differ significantly from those of Schwarzshcild holes only when J is very close to 
GM 2 /c. Our preliminary investigations of Kerr holes, in work now underway, confirms this. While 
frame dragging by the black hole may slightly increase the probability of detecting beams bent 
in the prograde sense, and decrease the probability of retrograde-bent beams, these effects largely 
cancel when considering the population as a whole. Thus the Schwarzschild approximation would 
appear to be justified for our purposes, in which it is the average properties that are of importance. 

We will rely heavily on results in Papers I and II for pulsar beam deflection around a Schwarzschild 
hole. In those papers (j)- m is the angle between the direction of pulsar emission and the direction 
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radially outward from the central SMBH at the emission event; the angle 4> ou t is the angle between 
that same radial direction and the direction in which the pulsar beam is moving when it is asymp- 
totically far from the SMBH, as sketched in Figure [H In the absence of the bending of the beam, 
the two angles (f>i n and ou t would be equal. The effect of curvature of the beam is encoded in the 
function F defined by 



^out 



F(<Pm,r ) 



(2) 



where ro is the distance of the emission point from the SMBH. The computational method for 
finding the F function is discussed in Papers I and II. A practical approximation for F, useful for 
the considerations of this paper, is presented in the Appendix. 
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pulsar spin 



axis 



orbital 
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Fig. 2. — The geometry of the orbit, pulsar spin, and pulsar beam. 



With the simplification to a nonrotating SMBH, we do not need to consider any angle between 
the orbit of the pulsar and the spin axis of the SMBH, The geometric parameters of interest are 
pictured in Fig. [2j The inclination of the pulsar spin axis with respect to the orbital plane is 
denoted A; the beam of pulsar emission is taken to have its center at angle «o from the spin axis, 
and to have width angular width 2Aa, so that the pulsar emission is confined between conical 
surfaces with opening angles ao — and «o + Aa, as shown in Fig. [2l 
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In Fig. [2j ro denotes the radial distance of the pulsar from the SMBH at the moment of 
emission of a beam. We do not assume circular oribits in our probability calculations except in the 
calculations of orbital times in Sec. HI 

We will assume that there is no favored alignment of the pulsar spin axis with the pulsar orbital 
plane, and will take A to be uniformly distributed over the sky. Within a fraction of a pc from 
the central SMBH, we expect there to be no significant alignment of the neutron star population 
with the disk of the Galaxy, so we will take the orientation of the orbital plane also to be randomly 
distributed over the sphere. 

We will assume the angle ato, between the pulsar beam and the spin axis to be randomly 
distributed over the half sphere, and we assume that for every choice of ao there is also a beam at 
7r — do- (In our probability estimates, we will avoid double counting in the case that ao + A a > vr/2 
and the conical regions of the two beams overlap.) 




n-urjn — ru- n diin- r 



0.1 0.2 0.3 

Pulse Duty Cycle (Fraction of Period) 



0.4 



Fig. 3. — A histogram of the FWHM of currently known pulsars in units of fraction of duty cycle. 

The angular width Act of the pulsar emission is a crucial parameter in the probability of obser- 
vation of deflected beams. Figure [3] shows a histogram, for currently known pulsars, of the FWHM 
of the pulse as a fraction of the pulse duty cycle. The mean of this distribution is 4.6%, so for prob- 
ability calculations in this paper we will use a duty cycle of 5%, and therefore a value of 9° for Aa. 

Of particular importance is the fact that a strongly bent beam will generally be reduced 
in intensity in comparison with a directly observed beam; thus our detection probabilities must 
account for the reduced brightness of the source. In terms of the bending function F of Paper I, 
the "amplification" factor (generally less than unity) for intensity is given by 

I_ = sin (f)- m 

I sin(F) (dF/d<f> in ) " U 
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As a step in understanding how much reduction can be allowed if a pulsar beam is to be observed, 
we use the radio luminosity at 1.4 GH^] for a ll pulsars for which thi s quan tity has been calculated 
in the Australia Telescope National Facility lATNF Pulsar Catalog! (|2010i ). We assume that the 
population of pulsars at the Galactic center has the same distribution of luminosities as these 
known pulsars. In that case, the distribution of pulsar flux densities S, observed at the Earth, 
would be that shown in Fig. [H 





£0 □ L 



1 
Log[S(mJy)] 



Fig. 4. — The distribution of pulsar L-band flux density S that would be seen on Earth if all 
cataloged pulsars were located at the center of the Galaxy. 

The minimum flux de tectable at a telescope can be estimated with the following equation from 



Lorimer &; Kramerl (J2004J) (hereafter LK): 



(SNR min )A)S sys / W 



Here W and P are respectively the pulsar pulse width and period. From our assumption that W/P 
is 5% we get yW/ (P — W) = 0.23. The parameter f3o is a correction factor for imperfections 
in data collection. Most current pulsar detection systems use a three-level correlarator with /3q of 
1.16 (see LK), and this is the value we shall use. For n p , the number of polarizations recorded and 
summed in the detection process, we will use n p = 2 because typically two polarizations are summed 
during pulsar detection scans. SNR m ; n is the minimum detectable signal-to- noise ratio required in 



The available data are mostly for this band; however, observations of beams from Sgr A* will likely need to be 
conducted in S band. Nonetheless, we expect the distribution of luminosities to be roughly the same in the two 
bands. 
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a search; we will take this to be 5. For the time pointed at the source, t Q b s , we will assume a 1 hour 
observing time. The bandwidth of the recorded data Af depends highly on the pulsar detection 
instruments used at a particular telescope. Bandwidths typically range from 100MHz to 800MHz. 
S sys is the system equivalent flux density which depends strongly on the collecting area of the 
telescope and the raw antenna sensitivity (see LK). The relevant characteristics of current existing 
radio telescopes and of possible future radio telescopes are detailed in Table CD 



Telescope Name 


^(Jy) 


A/(MHz) 


5 min (mJy) 


Parkes 


30 


340 


.022 


GBT 


10 


800 


.0048 


FAST 


1.5 


800 


.00072 


SKA 


.23 


800 


.00011 



Table 1: Theoretical S m i D values for two existing radio telescopes (Parkes, GBT) and two planned 
radio telescopes (FAST, SKA). The S sys values for the existing radio telescopes are measured system 
equivalent flux densities and the Af values are for current recording instruments. The S sys values 
for the two planned telescopes are theoretical values and the Af values were chosen to be the same 
as the best existing recording instruments. 



3. Probability calculations 

In this section we show how to calculate the probability that radiation from a single pulsar 
is detectable by a telescope on Earth, after having passed through the strong-field region of the 
black hole. This calculation naturally breaks down into two parts: determining what orientations 
of the pulsar and black hole produce strongly-bent beams, and determining where the Earth must 
be positioned relative to the system in order to detect those beams. 

In Paper II it was shown that for any relative position of pulsar, black hole, and Earth, there 
are a set of directions, called "keyholes," in which a photon could be emitted from the pulsar, 
pass around the black hole, and arrive at the Earth. These keyholes are typically within a few 
Schwarzschild radii of the black hole, so when the pulsar is far (many Schwarzschild radii) from 
the black hole, we can treat the keyhole as co-located with the black hole: that is, the pulsar beam 
must sweep across the black hole. The first part of the probability calculation is to determine for 
what fraction of the pulsar's orbit it is in a position to illuminate the black hole with its beam. 

We view the system from the perspective of the pulsar, so that the black hole traverses the sky 
of the pulsar along a great circle corresponding to the orbital plane. (This great circle in the pulsar 
sky does not imply that the pulsar-SMBH distance is constant.) Meanwhile, the pulsar spins about 
its rotation axis, and emits radiation in a cone offset from that axis: once per pulsar rotation the 
cone sweeps out an annulus in the sky of the pulsar. This is illustrated in Fig. [5j where ao, Aa 
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and A have the meaning described in Sec. [2] and pictured in Fig. [2j 

With Fig. [5jb) we introduce the angle (3, the total arc length (if any) over which the annulus 
intersects the orbital plane. When calculating (3, it is useful to focus on one hemisphere at a 
time and to label the two edges of the pulsar's radiation cone. We will define these two edges 
as aii = ao — Aa and «2 = «o + Aa. We break the calculation of (3 into three cases, with one 
case having two subcases. The first case is that in which the orbital plane of the system is never 
illuminated by the pulsar's radiation, 012 < A, as shown in Fig. [6l In this case, the value of (3 
is trivially zero. The second case is the case in which the orbital plane passes between the outer 
and inner edges of the annulus: cti < A < a<i- This case has two subcases. In the first subcase, 
that for «2 > 7r/2, the entire orbital plane is illuminated by the pulsar, as shown in Fig. [T^a), 
The value of (3 in this case is the entire range tt that lies in that hemisphere of the pulsar's sky. 
(Remember that we are assuming symmetric emission about the pulsar's rotational plane, so that 
the two beams together illuminate the full 2tt range of the orbital plane.) The second subcase, 
when «2 < 7r/2, has the outer edge of the annulus intersecting the orbital plane twice, as shown 
in Figure 0(b). In this case, (3 can be found from the spherical triangle version of Pythagoras's 
theorem, applied to the triangle with hypotenuse ai and sides A and (3/2: cos«2 = cosAcos/3/2, 
whence (3 = 2 arccos(cos 012/ cos A). The final case has both the outer and inner edges of the 
annulus crossing the orbital plane, cti > A, as illustrated in Figure El The calculation proceeds as 
in the previous case, but considers only the range of (3 between the triangles with hypotenuses 0.2 
and ol\\ (3 = 2 arccos(cos 02/ cos A) — 2 arccos(cos aij cos A). The cases and calculations for (3 are 
summarized in the following: 



Case 


Condition 


(3 Calculation 


1 

2a 
2b 
3 


«2 < A 

[ a 2 > vr/2 

OL\ < A < Q?2 < 

[ a 2 < vr/2 

A < cti 


(3 = 

(3 = TT 

/3 = 2arccos(^) 

/3 = 2arccos(^)-2arccos(^) 



The second part of the probability calculation is to determine how often the Earth will be in 
a position to detect strongly-bent beams from the system. This imposes geometric limitations to 
ensure that we are considering beams that are significantly deflected by the black hole, but that 
are not deflected so strongly that they are attenuated to a flux too low to be detected. These 
constraints turn out to place limits on the acceptable range of (pout- 

The determination of that range is shown in Fig. [9] for the case tq = 100M. That figure shows 
the dramatic amplification at </> ou t = TTj corresponding to strong lensing. For </> out slightly less than 
7T, the attenuation factor is unity, and the bending is not significant. There is bending for <j) oat < it, 
but the range of 4> ou t for which there is significant bending is small. Moreover, this range is even 
smaller than in Fig. [9] for the larger, more relevant values oiro/M. As a convenient approximation, 
therefore, W6 will consider "strong bending" only for (pout 

> tt. The figure shows that as 4> ou t 
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Fig. 5. — The geometry of a pulsar-black hole system as seen from the pulsar's reference frame. 
Panel (a) shows how the pulsar beam sweeps out an annular region in the pulsar's sky, which may 
intersect with the orbital plane of the pulsar-black hole system. Panel (b) shows the pulsar sky and 
labels the angles describing the orientation: A is the angle between the spin axis and orbital plane, 
ao the angle between the spin axis and the center of the illuminated annulus (or beam), and Aa is 
the half-width of the annulus (or beam). The zoom on the right shows the region in the pulsar's 
sky in which the Earth can lie if the photons reaching the Earth from the strong-field regime are 
to be bent by no more than an angle 5. 




Fig. 6. — Schematic of the pulsar sky for case 1, in which the orbital plane does not intersect the 
annulus. 
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Fig. 7. — Schematic of pulsar sky for case 2, in which the orbital plane passes between the outer 
and inner edges of the annulus. In panel (a), the orbital plane is entirely contained within the 
annulus. In panel (b), it crosses the outer edge of the annulus at two points symmetric about the 
meridian containing the spin axis. 




Fig. 8. — Schematic of pulsar sky for case 3, in which orbital plane intesects both the inner and 
outer edges of the annulus. 
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incresses beyond tt the attenuation becomes greater and greater. 

Our approach will be to specify a radius of emission ro and a minimum acceptable value of 
I /Iq. From a calculation like that shown in Fig. [9] we then find the value of the angle 5, the value 
of 4>out — 7T at which the attenuation is that of the minimum acceptable value of I /Iq. This value 
of 5 determines the range of directions in which the Earth must located if an Earth telescope is to 
detect the beam: the Earth must lie no more than an angle 5 from the pulsar-black hole axis. 

The corresponding region on the pulsar sky is illustrated in Fig. 0(b) . Since 5 is typically very 
small, we can express this area using a flat-space approximation: 2/3<5 + tt5 2 . The probability that 
observers on Earth can detect strongly-bent beams from a given pulsar is given by the size of this 
area over the angular area of one hemisphere (again, we assume symmetry across the rotation plane 
of the pulsar): 

P, = ,5) 

This of course assumes nonzero (3: if /3 = 0, the black hole is not illuminated by the pulsar, and 
there are no strongly-bent beams, so Pi = 0. Note that (3 is a function of the underlying parameters 
A, an d Aa, and note that 5 is a function of the underlying parameters vq/M and I m i n /Io- 




Fig. 9. — The intensity ratio I /Iq as a function of </> ou t) for ro = 100M. The shaded regions show 
the ranges of </> ou t; that are not of interest either because the pulsar beam is too weakly deflected 
or because it is too strongly attenuated. 
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4. Results 

Now that we have shown how to calculate the probability of Earth detection for a particular 
pulsar, we will describe how we estimated the number of pulsars that would be detected given 
assumptions about the distributions of pulsar characteristics. We ran Monte Carlo simulations 
that selected A from a uniform distribution of cos A, so that the direction of the pulsar spin axis 
was uniformly distributed over the sky; similary ao was chosen from a uniform distribution of 
cosao- Then a value for the pulsar's flux, S, was chosen from the distribution shown in Fig. SJ 
Lastly, a value of tq/M was chosen from the distribution in Eq. (pQ); we cut this distribution off 
at tq/M = 400,000 since pulsars beyond that point contribute little to the total probability (see 
below). 

The simulations took the pulsar parameters (ro, S) chosen by the Monte Carlo method. From 
these a determination was made of the minimum value of I /I$ (equivalently S m i n /S) that can be 
detected for those pulsar parameters. From 1/1$ and tq, the value of 5 was determined. The value 
of (3 was detemined through the calculations described in the previous section applied to the values 
of A and ao chosen by the Monte Carlo method. The justification for using Aa = 9° has been 
explained in Sec. [2J Pi was then calculated with Eq. [SJ This gives us the probability that a single 
pulsar's beam will explore the black hole's strong gravitational field and still be detectable once it 
reaches the Earth. Monte Carlo simulations were repeated to insure that results were consistent to 
better than 1%. We then multiplied our result by the total number of pulsars out to ro = 400, 000, 
according to the Pfahl and Loeb distribution of Eq. ([TJ . The result is the total number of pulsars 
(-ftot) that will, at some point in their orbit about the central SMBH, both illuminate the strong 
gravitational field and be detectable at Earth. Table [2] gives this number for the four telescopes 
considered in Sec. [2j 



Telescope Name 


'S'sys(Jy) 


-ftot 


Parkes 


.022 


8.75776 


GBT 


.0048 


15.8075 


FAST 


.00072 


111.677 


SKA 


.00011 


207.021 



Table 2: S m i Q and Ptot values for two existing radio telescopes (Parkes, GBT) and two planned 
radio telescopes (FAST, SKA). 

The number of pulsars that are observable at some point in their orbit is not directly relevant 
if the orbital time is much larger than the duration of an observing program. For that reason we 
introduce a more useful number, the "observability," P bs 5 to describe the probability of observing 
a given pulsar in a limited-time observing program. In order to calculate this number, we ran 
the Monte Carlo simulations with a specified observational program duration (L). Once Pi was 
calculated, we then compared the pulsar's orbital period (T) to L. If T was less than or equal to 
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L, then P t, s was taken to be P\. If T was greater than L, then P\ was replaced by 

Pi x L 



obs 



(6) 



and the result was multiplied by the total number of pulsars. Figure [10] shows the results of the 
Monte Carlo simulations for observational program durations ranging from one year to seven years. 
Figure [11] shows the variation in the number of pulsars detected, for S m { n = 0.022 mJY, as the 




1 2 3 4 5i 

Observational Program Length(Y) 



Fig. 10. — Expected number of pulsars observed as a function of observational program duration. 
Curves are labeled with the values of S m i n , in mJy, assumed in the computations. 



cutoff radius is changed, and justifies our use of the cutoff at r /M = 4 x 10 5 . 




2e+05 4e+05 6e+05 8e+05 le+06 
Cutoff value of r Q /M 



Fig. 11. — The number of pulsars expected, for S m i n = 0.022 mJy, as a function of the cutoff tq/M 
for the distribution in Eq. ([!]). Each curve is labeled with the number of years assumed for the 
search. 
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5. Observing Programs and Strategies 

A natural first question about observing strongly deflected beams is "how will we know that 
they are strongly deflected?" The answer starts with the fact that the angle through which the 
beam is "strongly" deflected is not large. For our paradigmatic case, tq = 10 4 M, the bending is 
approximately 0.036 rad. For larger tq the deflection, for a given S^m^, is even smaller. 

Since the "strong" deflection is small, we will receive a deflected beam only when the emitting 
pulsar, the SMBH and the Earth are almost on a straight line. Since pulsar beam widths are large 
compared to the deflection, this means that if the Earth receives the deflected beam, it will also 
receive the direct beam. The geometry of the direct and deflected beams is shown in Fig. [T2l where 
we see that the angle, at reception, of the direct and deflected beams is not of order tq divided by 
the Earth-SMBH distance (8 kpc), but rather is of order of that number multiplied by the deflection 
angle. The result, an angle of order 10 -8 radians, is less than the resolution of radio telescopes. 
We conclude that any monitoring of the innermost region of Sgr A* for a deflected beam, will also 
monitor for a direct beam. 

The criterion for the detection of a strongly deflected beam will be that it is one of a pair of 
pulses detected with very similar pulse periods, differing only due to a phase modulation caused by 
variations in the propagation times along the two paths: the two sets of pulses will differ in pulse 
period by a fractional amount of order the pulsar velocity divided by c. (See Papers I and II for 
more detail on phase effects and intensity effects of deflection.) If a deflected beam is detected we 
therefore assume that it will be relatively simple to identify it as deflected. 



deflection deflected beam 




pulsar 

Fig. 12. — The geometry of the observation of direct and deflected beams. 

If the Pfahl and Loeb distribution of Eq. ([I]) is approximately valid, the results of the previous 
section, especially Fig. [TUJ suggest that there is a very good chance of observing a strongly deflected 
pulsar beam with an observing program of 3-5 years, using existing telescopes. Monitoring of 
Sgr A*, of course, will face the problem of dispersion by the plasma density in the Galactic center. 
Assuming dispersion measures in excess of 1000pc/cm 3 , searches would have to be conducted in S- 
band (2-4 GHz) or higher frequencies in order to resolve pulses, rather than the more conventional 
L-band. 

We note that the pulsar fluxes and current telescope sensitivities were computed for the L 
band, since that is where most pulsar luminosities have been measured. Relative sensitivities in 
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the S band may be somewhat lower. However, we note that the baseline value of S m in = 0.022 mJy 
for the Parkes telescope is very conservative in view of the sensitivities that will soon be available. 
See Table [IJ and the corresponding estimates of detectable pulsars in Fig. [TOj Even with the 
current estimated senstitivity of Parkes there is a good chance of detecting strongly-bent pulsar 
beams in a 2-3 year observing program; with more sensitive telescopes such as GBT or FAST 
(under development), the number could rise to several pulsars. Longer observing times increase the 
number of detectable pulsars by allowing long-period pulsars more time to come into alignment. 

An appropriate multiyear observing program can be carried out "in background" at a telescope. 
Observations must be made sufficiently frequently not to miss the relatively short epoch during 
which the pulsar/SMBH/Earth alignment leads to a strongly deflected beam meeting the I/Iq 
criterion. That epoch is of order of the orbital time multiplied by the ratio 5/(2ir). For our 
prototypical choice ro = 10 4 M, this epoch is on the order of a week. Searches for deflected beams 
in Sgr A*, therefore, would have to be made every other day. The observing session would be of 
a duration of that used for any other pulsar search, on the order of 1 hour. (Longer observations 
might be considered, in view of the effect shown in Eq. (Jl|) of t b s on 5 m i n , and hence on the I/Iq 
cutoff. 

An important question to ask is why we do not yet have evidence of the assumed large popu- 
lation of pulsars in Sgr A* . This can be explained by the lack of any concerted effort to search the 
Galactic center for pulsars, and the need for S band searches with large dispersion corrections. In 
addition the Galactic Center is out of view for the Arecibo telescope, and in a portion of the sky 
with reduced sensitivity for GBT. Currently, Parkes is the largest telescope with a clear view of 
the Galactic Center; the FAST telescope will substantially increase our sensitivity to the Galactic 
Center. 

It should be noted that, at least to some extent, monitoring the Galactic center for strongly 
deflected beam would constitute a more general search for pulsars in that region. The possibility 
of detecting deflected beams and making measurements of the parameters of the SMBH provides 
added scientific motivation for such a survey. 

6. Conclusions 

Our estimates suggest that a multi-year program that monitors Sgr A* with radio observations 
for one hour every other day has a reasonable probability of detecting pulsar beams that have 
been strongly deflected by our Galaxy's SMBH. With instruments coming in the near future, 
in particular, FAST and SKA, the probability should become high enough so that a three year 
observational program either detects a strongly deflected beam, or that the failure to make such a 
detection puts useful limits on the density of pulsars in the Galactic Center. 

Our estimates in this paper constitute a first step in the study of probabilities of detection of 
a strongly deflected beam. The intention was to establish whether the probabilities are so small 
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that ob servations are out of the question, or so large that current observations rule out models, like 



that of lPfahl fc Loebl (|2004l ). with a significant density of pulsars in Sgr A*. The estimates in this 
paper establish neither extreme: a concerted observing program with the best current telescopes 
would optimistically detect only a few such pulsars. This provides motivation for such a program, 
and also for further study of the problem of pulsar beam deflection by SMBHs. 

Such an improved study would have to include effects of spin of the SMBH, and of eccentricity 
of orbits. While our approach of using averages and simple assumptions was appropriate to the 
purpose of this paper, effects due to SMBH spin, and high eccentricity, could increase the parameter 
space of pulsar configurations whose beams can reach the Earth. Such work is now underway. 

The most exciting result of these preliminary estimates is their indication that we are poten- 
tially on the verge of detecting a new phenomenon: pulsar beams that have passed through the 
strong field region of the SMBH at the center of our Galaxy, beams that can bring us information 
about the properties of the SMBH and its surrounding spacetime might be inaccessible in any other 
way. 
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PHY0554367 and HRD0734800. We also thank the Center for Gravitational Wave Astronomy at 
the University of Texas at Brownsville. 



7. Appendix 

We are primarily interested in values of <p\ n that are only slightly smaller than it. In this case 
the photon path starting at some very large radius will penetrate to small radii, and almost all the 
bending will take place at small radii. We can then find the bending for emission from the very 
large radius, say Rq = 10, 000 M, by considering the bending only interior to a smaller large radius, 
say ro = 100 M. In effect, we are considering the large radius region from R Q to r G to be flat space. 
From the the curve F((j)[ n ; r) for ro = 100 M, therefore, we can infer the curve for Rq = 10, 000 M, 
and for all larger radii (provided, of course, that (f)- m is near tt so that both R Q and r D are much 
larger than the radii at which the bending occurs). 

To find the 4>^^ for a photon emitted at a very large Ro, we choose the smaller ro along the 
future path of the photon to be a radius for which we know the curve = ^(^ui > r o)- Here 
superscripts (ro) and {Ro) distinguish the angles associated with the two radii. The Euclidean 
geometry relating </>;f and (/>| n ^ ^ is given by the law of sines to be 

^^sm-^iRo/rojsm^) , (7) 

and we choose the branch of sin -1 so that <^ o) > n/2. 

We next notice that </> ou t for the photon starting at ro is less than out for that same photon 
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Fig. 13. — The flat spacetime propagation from large radius to very large radius. Angles <p- m for 
the large, and the very large emission radius are distinguished with the superscripts (ro) and (Ro), 
respectively. 



world line considered to start at Rq, in flat space, according to 

JRo) _ An) 



J out 



'out 



+ 7- 



With 7 evaluated in terms of the ingoing angles, this becomes 



6 (Ro) _ ,(ro) ,(Ro) 
''out — Vout ' V^in 



(8) 



(9) 



From a combination of Eqs. (J7|) and ([9]), written in terms of the function F((j>\ r), the full expression 
can be given as: 



F(<j>;Ro) 



+ F 



sin 1 [^ a sini 

ro 



; r 



sul _1 [ sin< 

ro 



(10) 



Thus, knowing the F function for any (sufficiently large) ro, we can evaluate it for any larger Rq, 
and thus determine the maximum deflection angle 5 via Eq. ([2D . Figure [J3] shows our reference 
bending function for a radius ro = 100M. 
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